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I. INTRODUCTION 

Entanglement is a unique feature of quantum systems 
that play a key role in quantum information processing. 
The early study of entanglement is only focused on the 
foundations of the quantum mechanics [l| . Recently due 
to its potential applications in quantum communications, 
quantum cryptography, quantum computer and quantum 
information [2j, entanglement has been studied exten- 
sively. One of the most important progress is the quanti- 
tative measures of entanglement for mixed state by using 
the entanglement of formation d H[ . For the special case 
of two spin-1/2 systems, the entanglement of formation 
is given by the concurrence [H, Q. Newly, considerable 
interest has been devoted to entanglement of quantum 
spin system 0, H|, identical particles [E EE], fractional 
quantum Hall effect [ll[, and spins of a noninteracting 
electron gas fl2l ]. 

On the other hand, since Anderson published his fa- 
mous paper [1 31 ] about disorder induced localization, ex- 
tensive investigations have focused on the metal-insulator 
transition (MIT). For one-dimensional (ID) Anderson 
model, it is well known ,14:] that all eigenstates are lo- 
calized and there is no mobility edge separating local- 
ized and extended states. However, the specific extended 
states and/or mobility edges have been found in sev- 
eral ID determined and random model with short-range 
and long-range correlation [IE EE El El El EE El- 
The well-studied examples of the potentials are a slowly 
varying potential [H EE El , random-dimer potential 
[HI], long-range correlated disordered potential fis)l and 
Anderson model with long-range hopping [20l |2l| . In 
these models, electronic localized behaviors are studied 
by judging Thouless exponent (or Lyapunov coefficient), 
participation ratio or dynamics of wave function. 

Recently mode entanglement of spinless electrons shar- 
ing in one-particle states in ID model has been investi- 
gated in Refs. [H EH- Using the ordinary Harper and 



the kicked Harper model, Lakshminarayan and Subrah- 
manyam found that entanglement can reflect MIT. Simi- 
lar behavior is also found for the ground state of an elec- 
tron in ID Frenkel-Kontorova (FK) potential [H[. There 
are many more complex ID potentials, e.g., potentials 
used in Refs. 0, EE El El, EE HE El , which have been 
widely used to study the MIT and exhibit more complex 
localization behaviors than that of the Harper model. 
Therefore, it is interesting to study the mode entangle- 
ment of an electron in these more complex ID potentials. 

The paper is organized as follows. In next section the 
formalism of entanglement and concurrence is described. 
In Sec. Mil the numerical results for four kinds of models 
are presented. Section IIVI is devoted to conclusion. 



II. FORMALISM 

In the second-quantized picture, the Hamiltonian for 
electrons moving in ID determined and/or random po- 
tential can be written as follow: 
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H = -t^2 ( C n C n+l + 4+l C «) + V ^ C t°n, (1) 

n—1 n—1 

where t is a nearest-neighbor hopping integral, c+ (c n ) is 
the creation(annihilation) operator of nth site, and V n is 
the one-site potential. In our numerical studies, we take 
t = 1 without loss of generality. The site occupation 
basis is 
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where rt,; = 0, 1, and |0) is the vacuum. Note that there 
is an isomorphism between these states and the states of 
N qubits[22|]. For an electron, J2 
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is 



|0,. 



= 1 . If we write 
, 0) , the general state of an electron 
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where ty n is amplitude of wave at nth site. 
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From Eqs. (1), (2) and (3), we obtain the eigenequa- 
tion 



(4) 



where E is the eigenenergy. For an eigenstate f3 with 
eigenenergy Ep, the concurrence between sites ( or 
qubits) i and j is given [22j as 



C^2 



« j 



(5) 



States that have a large minimum pairwise concurrence 
can be said share entanglement better. Specifically, when 

tyf = 1/VJV, the state becomes the so-called W state 

24j and the concurrence is given by 2/N. There can 
not be states whose minimum pairwise concurrence ex- 
ceeds 2/N. As a gross but useful measure of entangle- 
ment sharing, Lakshminarayan and Subrahmanyam pro- 
pose and study the pairwise concurrence averagely in a 
given state. For the given eigenstate, 



(6) 



i<j 



where d — N(N — l)/2. From the definition of (6), we 
can see that < C@ > has connections to measures of 
localization of the eigenstates. As a further gross measure 
they also average over all the eigenstates f3 , i.e. , 



(7) 



where M is the number of all the states. In the following 
concurrence (C) and (C^^ are measured for four kinds 
of ID determined and random potentials . 



III. NUMERICAL RESULTS 
A. Slowly varying potential 

In the slowly varying potential model [HI, [l7| , the 
on-site potential is given by V n — Xcos(nan v ), here A, a 
and < v < 1 are positive numbers which completely de- 
fine the tight-binding problem. If v = 1, the model is well 
known Harper model. For Harper model, the eigenstates 
are either all extended or all localized states depending 
on whether A is smaller or larger than 2.0 [l7|. When 
A = 2.0, all the states are critical. The mode entangle- 
ment of Harper model has been investigated in Ref. |22| 
and it was found a sharp transition in concurrence at 
A = 2.0, which corresponds to MIT. 

For < v < 1, it is well known that there are two 
mobility edges at E c — ±(2.0 — A) provided that A < 2.0. 
It was found extended states in the middle of the band 
( \E\ < 2.0 — A ) and localized states at the band edge 
( 2.0 - A < \E\ < 2.0 + A ). For A > 2.0, all states are 
found to be localized. 
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FIG. 1: Average concurrence N (C 13 ) of the individual state 
in the slowly varying potential model as functions of energy, 
with 7ra = 0.2 and v = 0.7, (a) A = 0.4, (b) A = 2.0 and (c) 
A = 4.0, respectively, (d) Average concurrence JV (C) over all 
states varying with A. 
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Concurrence N (C 3 ) are plotted in Fig. 1(a), (b), (c) 
for A = 0.4,2.0 and 4.0 respectively. From Fig. 1(a), (b), 
it's clearly shown that there are sharp transitions in con- 
currence at the mobility edges E c = ±(2.0 — A). The 
transition becomes sharper when N increasing. At ex- 
tended states, N (C^) w 1.6. At the localized states, 
N (C^) is decreasing from mobility edges to band top 
(or bottom). In Fig.l (c), all states are localized and 
N (C' 3 ) are small for all states, which is smaller than 
1.6. Comparing with Fig. 1(a), (b) and (c) in Ref . [l7j . 
we can found that the longer the localized length is, the 
larger concurrence is. 

The concurrence N (C) averaged over all states is plot- 
ted in Fig. 1(d). There is a transition at A = 2.0. The 
transition becomes sharper when N increasing. Obvi- 
ously there is dramatically transition as N — > oo at 
A = 2.0. 



B. Random-dimer potential 

Another interesting ID model, which has localized and 
extended states, is random-dimer model [l8l |. In this 
model, the site energies V a and Vb are assigned at random 
to 2nth (n=integer) site with probability q and l — q, and 
V2n+i = V^n- By solving the time-depended Schrodinger 
equation, Dunlap et al. [T§ | found that the mean-square 
displacement at long times is shown to grow in time as 
i 3 / 2 provided -2 < V a - V b < 2, diffusion occurs if 
V a — Vb = ±2 and localization otherwise. These mean 
that extended states exist when —2 < V a — Vb < 2, and 
there are only localized states when \V a — Vb\ > 2. 

Concurrence iV^C' 3 ) are plotted in Fig. 2(a), (b), (c) 
for V a - V b = 1.0, 2.0 and 2.5 respectively. The V b is 
taken as 1.0 without loss of generality and q is equal to 
0.5 corresponding to most random situation. The results 
of Fig. 2 are obtained for average of 200 samples [25|. The 
averages with more samples give same results. There are 
two bumps in concurrence when V a — Vj, = 1.0. There 
is no obvious mobility edge in this model, so there is no 
sharp transition in concurrence as that in the slowly vary- 
ing potential model. There are two jumps in concurrence 
when V a — V b — 2.0. For V a — Vb = 2.5, concurrences are 
small for all states. 

Average concurrence N (C) as functions of V a — Vb is 
plotted in Fig. 2(d). There is a jump in the concurrence 
at V a — Vb = 2.0, which is in accordance with the critical 
value of V a — V> obtained by dynamical method. 



C. Long-range correlated disordered potential 

Recently, another kinds of disordered potential is stud- 
ied extensively. The potential is self-affine Gaussian po- 
tential with Hurst exponent < H < 1, such that 
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FIG. 2: Average concurrence iV (C 3 ) of the individual state in 
the random-dimer model as functions of energy, (a) V a — Vb = 
1.0, (b) V a -Vb = 2.0 and (c) V a -V b = 2.5, respectively, (d) 
Average concurrence JV (C) over all states in the dimer model 
as functions of V a — Vb. 



tential can be generated by fractional Brownian motion. 
To generate the trace of a fractional Brownian motion, 
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an approach based on the use of discrete Fourier trans- 
forms to construct such long-range correlated sequences 
can be applied. The on-site energies can be given by the 
relation tl9|: 



N/2 

fc=l 



2tt 



-a) 1/2 



cos( 
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where TV is the number of sites and tfk are TV/2 indepen- 
dent random phases uniformly distributed in the interval 
[0, 2ir]. The sequence usually has an approximate power- 
law spectral density of the form S(k) oc l/k a , where 
S(k) is the Fourier transform of the two-point correla- 
tion < ViVj >. Here a = 2H + 1. Same as in Ref. [H|, 



we normalized the energy sequence to have < V n > 
and AV -- 





(Vn) ~ (Vn) 2 = 1. 
It is shown [l9( that when a < 2.0, all states are lo- 
calized; when a > 2.0, localized states occur in the edges 
of the band and the extended states in the middle of the 
band, separated by mobility edge. 

The relations between concurrence TV(C^) and 
eigenenergy E are plotted in Fig. 3 (a), (b) and (c) for 
a = 5.0,2.0 and 1.0 respectively. Here 200 samples of 
random iff. are averag ed. For a = 5.0, N (C 13 ) « 1.6 
in the middle states and N (C^j is gradually decreas- 
ing from center to band top (or band bottom). The 
eigenenergies region of extended states is decreasing as 
a becomes smaller [l9j]. When a = 2.0, the concurrence 
is near 1.6 only at band center, which is shown clearly 
in inset of Fig. 3(b). When a — 1.0, concurrence for all 
states is small, which corresponds to localized states. 

Concurrence averaged over all states is plotted in 
Fig. 3(d). There is an inflexion when a is near 2.0 in the 
figure. Same as in model with slowly varying potential, 
the bigger TV is, the transition is sharper. 



D. Random potential with long-range hopping 

The random potential with long-range hopping has re- 
ceived considerable attention recently. The Hamiltonian 
of such ID tight-binding model is expressed as 



H 



^e„ \n) (n\ + 

^ ^ J inn 

n n^m 



(9) 



where e n is energy level at nth site , uniformly distributed 



in interval [-W/2, W/2] and J mn = J/\m- 



(Jn 



0) is the long-range hopping amplitude. We will adopt J 
as energy units without loss of generality. The periodic 
boundary condition is applied. 

By using a supersymmctric method combined with a 
renormalization group analysis, Rodriguez et al. (2p| 
have shown the existence of extended states for energies 
within a range near the band top in one and two dimen- 
sional Anderson models. They found that MIT occurs 
only within the range of d < /i < 3d/2 in the thermo- 
dynamic limit, no matter how large the value of W is, 
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FIG. 3: Average concurrence iV (C 3 ) of the individual state in 
the long-range correlated disordered potential model as func- 
tions of energy for (a) a = 5.0, (b) a = 2.0 and (c) a = 1.0, 
respectively. The insets show concurrence for a typical dis- 
order realization at N=800. (d) Average concurrence N (C) 
over all states as functions of a. In all figures, from up to low, 
TV = 200, 400 and 800. 



5 



here d is the geometric dimensionality of the system. Us- 
ing finite size scaling analysis combined with the transfer 
matrix method, Xiong and Zhang [2l| found that there 
exists MIT at critical value \i c for some W. 

The concurrence N (C^^ as functions of energy are 
plotted in Fig. 4 (a), (b) and (c) for fi = 1.1, 1.5 and 
1.7 respectively. Here W — 5 is taken as example and 
the phenomena are similar for other W. The results 
are obtained for average of 200, 100, 50 samples [25[ for 
N = 200,400,800 respectively. When fi = 1.1 and 1.5, 
concurrence N (C 3 ) is large for states near the band top, 
which means there exist extended states. It is small for 
the states near the band bottom. This is quite different 
from that of above three models. For /i = 1.7, the con- 
currence is small for all states, which means all states are 
localized. 

The concurrence averaged over all states is plotted in 
Fig. 4(d). From the inset, we can get the inflexion for /i 
near 1.70, which is consistence with the upper limit for 
critical value /x c obtained in Ref.[2lJ. 

IV. CONCLUSION 

Using the measure of concurrence, mode entanglement 
sharing in one-particle states in four kinds of models 
is studied numerically. Concurrence N (C' 3 ) at a given 
state and N (C) averaged over all states are investigated. 
For (C' 5 ), the concurrence is large in extended states 
and small in localized states. There is a sharp transi- 
tion in the concurrence at mobility edge. N (C 3 ) gives 
the information about the localization behavior of the 
given eigenstate f3. From the curves of the N (C) vs 
the parameter of the models, which is A, V a — V&, a, 
and /i for slowly varying potential, random-dimer po- 
tential, long-range correlated disordered potential, and 
long-range hopping random potential, respectively, we 
can found clearly that there is an inflexion (or jump) at 
a critical parameter value, which is in accordance with 
that obtained by other methods. When parameter value 
is greater (smaller) than the critical parameter value, the 
system has only localized eigenstates, while when param- 
eter value is smaller (greater) than the critical parame- 
ter value, the system has both localized and delocalized 
states, which is different from that of one-dimensional 
Harper model. The inflexion or transition point in the 
curve of N (C) versus parameter of systems corresponds 
to the disappear of delocalized states. Therefore mode 
entanglement can be a new index to reflect MIT. 
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I. INTRODUCTION 

Entanglement is a unique feature of quantum systems 
that play a key role in quantum information processing. 
The early study of entanglement is only focused on the 
foundations of the quantum mechanics [l| . Recently due 
to its potential applications in quantum communications, 
quantum cryptography, quantum computer and quantum 
information [2j, entanglement has been studied exten- 
sively. One of the most important progress is the quanti- 
tative measures of entanglement for mixed state by using 
the entanglement of formation d H[ . For the special case 
of two spin-1/2 systems, the entanglement of formation 
is given by the concurrence [H, Q. Newly, considerable 
interest has been devoted to entanglement of quantum 
spin system 0, H|, identical particles [E EE], fractional 
quantum Hall effect [ll[, and spins of a noninteracting 
electron gas fl2l ]. 

On the other hand, since Anderson published his fa- 
mous paper [1 31 ] about disorder induced localization, ex- 
tensive investigations have focused on the metal-insulator 
transition (MIT). For one-dimensional (ID) Anderson 
model, it is well known ,14:] that all eigenstates are lo- 
calized and there is no mobility edge separating local- 
ized and extended states. However, the specific extended 
states and/or mobility edges have been found in sev- 
eral ID determined and random model with short-range 
and long-range correlation [IE EE El El El EE El- 
The well-studied examples of the potentials are a slowly 
varying potential [H EE El , random-dimer potential 
[HI], long-range correlated disordered potential fis)l and 
Anderson model with long-range hopping [20l |2l| . In 
these models, electronic localized behaviors are studied 
by judging Thouless exponent (or Lyapunov coefficient), 
participation ratio or dynamics of wave function. 

Recently mode entanglement of spinless electrons shar- 
ing in one-particle states in ID model has been investi- 
gated in Refs. [H EH- Using the ordinary Harper and 



the kicked Harper model, Lakshminarayan and Subrah- 
manyam found that entanglement can reflect MIT. Simi- 
lar behavior is also found for the ground state of an elec- 
tron in ID Frenkel-Kontorova (FK) potential [H[. There 
are many more complex ID potentials, e.g., potentials 
used in Refs. 0, EE El El, EE HE El , which have been 
widely used to study the MIT and exhibit more complex 
localization behaviors than that of the Harper model. 
Therefore, it is interesting to study the mode entangle- 
ment of an electron in these more complex ID potentials. 

The paper is organized as follows. In next section the 
formalism of entanglement and concurrence is described. 
In Sec. Mil the numerical results for four kinds of models 
are presented. Section IIVI is devoted to conclusion. 



II. FORMALISM 

In the second-quantized picture, the Hamiltonian for 
electrons moving in ID determined and/or random po- 
tential can be written as follow: 
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H = -t^2 ( C n C n+l + 4+l C «) + V ^ C t°n, (1) 

n—1 n—1 

where t is a nearest-neighbor hopping integral, c+ (c n ) is 
the creation(annihilation) operator of nth site, and V n is 
the one-site potential. In our numerical studies, we take 
t = 1 without loss of generality. The site occupation 
basis is 
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where rt,; = 0, 1, and |0) is the vacuum. Note that there 
is an isomorphism between these states and the states of 
N qubits[22|]. For an electron, J2 
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= 1 . If we write 
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where ty n is amplitude of wave at nth site. 
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From Eqs. (1), (2) and (3), we obtain the eigenequa- 
tion 



(4) 



where E is the eigenenergy. For an eigenstate f3 with 
eigenenergy E, the concurrence between sites ( or qubits) 
i and j is given [22j as 



C^2 



(5) 



States that have a large minimum pairwise concurrence 
can be said share entanglement better. Specifically, when 

tyf = 1/VJV, the state becomes the so-called W state 

24| and the concurrence is given by 2/N. There can 
not be states whose minimum pairwise concurrence ex- 
ceeds 2/N. As a gross but useful measure of entangle- 
ment sharing, Lakshminarayan and Subrahmanyam pro- 
pose and study the pairwise concurrence averagely in a 
given state. For the given eigenstate, 



(6) 



i<j 



where d — N(N — l)/2. From the definition of (6), we 
can see that < C@ > has connections to measures of 
localization of the eigenstates. As a further gross measure 
they also average over all the eigenstates f3 , i.e. , 



(7) 



where M is the number of all the states. In the following 
concurrence (C) and (C^^ are measured for four kinds 
of ID determined and random potentials . 



III. NUMERICAL RESULTS 
A. Slowly varying potential 

In the slowly varying potential model [HI, [l7| , the 
on-site potential is given by V n — Xcos(nan v ), here A, a 
and < v < 1 are positive numbers which completely de- 
fine the tight-binding problem. If v = 1, the model is well 
known Harper model. For Harper model, the eigenstates 
are either all extended or all localized states depending 
on whether A is smaller or larger than 2.0 [l7|. When 
A = 2.0, all the states are critical. The mode entangle- 
ment of Harper model has been investigated in Ref. |22| 
and it was found a sharp transition in concurrence at 
A = 2.0, which corresponds to MIT. 

For < v < 1, it is well known that there are two 
mobility edges at E c — ±(2.0 — A) provided that A < 2.0. 
It was found extended states in the middle of the band 
( \E\ < 2.0 — A ) and localized states at the band edge 
( 2.0 - A < \E\ < 2.0 + A ). For A > 2.0, all states are 
found to be localized. 
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FIG. 1: Average concurrence N (C 13 ) of the individual state 
in the slowly varying potential model as functions of energy, 
with 7ra = 0.2 and v = 0.7, (a) A = 0.4, (b) A = 2.0 and (c) 
A = 4.0, respectively, (d) Average concurrence JV (C) over all 
states varying with A. 
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Concurrence N (C 3 ) are plotted in Fig. 1(a), (b), (c) 
for A = 0.4,2.0 and 4.0 respectively. From Fig. 1(a), (b), 
it's clearly shown that there are sharp transitions in con- 
currence at the mobility edges E c = ±(2.0 — A). The 
transition becomes sharper when N increasing. At ex- 
tended states, N (C^) w 1.6. At the localized states, 
N (C^) is decreasing from mobility edges to band top 
(or bottom). In Fig.l (c), all states are localized and 
N (C' 3 ) are small for all states, which is smaller than 
1.6. Comparing with Fig. 1(a), (b) and (c) in Ref . [l7j . 
we can found that the longer the localized length is, the 
larger concurrence is. 

The concurrence N (C) averaged over all states is plot- 
ted in Fig. 1(d). There is a transition at A = 2.0. The 
transition becomes sharper when N increasing. Obvi- 
ously there is dramatically transition as N — > oo at 
A = 2.0. 



B. Random-dimer potential 

Another interesting ID model, which has localized and 
extended states, is random-dimer model [l8l |. In this 
model, the site energies V a and Vb are assigned at random 
to 2nth (n=integer) site with probability q and l — q, and 
V2n+i = V^n- By solving the time-depended Schrodinger 
equation, Dunlap et al. [T§ | found that the mean-square 
displacement at long times is shown to grow in time as 
i 3 / 2 provided -2 < V a - V b < 2, diffusion occurs if 
V a — Vb = ±2 and localization otherwise. These mean 
that extended states exist when —2 < V a — Vb < 2, and 
there are only localized states when \V a — Vb\ > 2. 

Concurrence iV^C' 3 ) are plotted in Fig. 2(a), (b), (c) 
for V a - V b = 1.0, 2.0 and 2.5 respectively. The V b is 
taken as 1.0 without loss of generality and q is equal to 
0.5 corresponding to most random situation. The results 
of Fig. 2 are obtained for average of 200 samples [25|. The 
averages with more samples give same results. There are 
two bumps in concurrence when V a — Vj, = 1.0. There 
is no obvious mobility edge in this model, so there is no 
sharp transition in concurrence as that in the slowly vary- 
ing potential model. There are two jumps in concurrence 
when V a — V b — 2.0. For V a — Vb = 2.5, concurrences are 
small for all states. 

Average concurrence N (C) as functions of V a — Vb is 
plotted in Fig. 2(d). There is a jump in the concurrence 
at V a — Vb = 2.0, which is in accordance with the critical 
value of V a — V> obtained by dynamical method. 



C. Long-range correlated disordered potential 

Recently, another kinds of disordered potential is stud- 
ied extensively. The potential is self-affine Gaussian po- 
tential with Hurst exponent < H < 1, such that 



< {Vn, 



V n 



> = 



:2H 



Such sequence of po- 
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FIG. 2: Average concurrence iV (C 3 ) of the individual state in 
the random-dimer model as functions of energy, (a) V a — Vb = 
1.0, (b) V a -Vb = 2.0 and (c) V a -V b = 2.5, respectively, (d) 
Average concurrence JV (C) over all states in the dimer model 
as functions of V a — Vb. 



tential can be generated by fractional Brownian motion. 
To generate the trace of a fractional Brownian motion, 
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an approach based on the use of discrete Fourier trans- 
forms to construct such long-range correlated sequences 
can be applied. The on-site energies can be given by the 
relation tl9|: 



N/2 

fc=l 



2tt 



-a) 1/2 



cos( 



2nik 
N 



(8) 



where TV is the number of sites and tfk are TV/2 indepen- 
dent random phases uniformly distributed in the interval 
[0, 2ir]. The sequence usually has an approximate power- 
law spectral density of the form S(k) oc l/k a , where 
S(k) is the Fourier transform of the two-point correla- 
tion < ViVj >. Here a = 2H + 1. Same as in Ref. [H|, 



we normalized the energy sequence to have < V n > 
and AV -- 





(Vn) ~ (Vn) 2 = 1. 
It is shown [l9( that when a < 2.0, all states are lo- 
calized; when a > 2.0, localized states occur in the edges 
of the band and the extended states in the middle of the 
band, separated by mobility edge. 

The relations between concurrence TV(C^) and 
eigenenergy E are plotted in Fig. 3 (a), (b) and (c) for 
a = 5.0,2.0 and 1.0 respectively. Here 200 samples of 
random iff. are averag ed. For a = 5.0, N (C 13 ) « 1.6 
in the middle states and N (C^j is gradually decreas- 
ing from center to band top (or band bottom). The 
eigenenergies region of extended states is decreasing as 
a becomes smaller [l9j]. When a = 2.0, the concurrence 
is near 1.6 only at band center, which is shown clearly 
in inset of Fig. 3(b). When a — 1.0, concurrence for all 
states is small, which corresponds to localized states. 

Concurrence averaged over all states is plotted in 
Fig. 3(d). There is an inflexion when a is near 2.0 in the 
figure. Same as in model with slowly varying potential, 
the bigger TV is, the transition is sharper. 



D. Random potential with long-range hopping 

The random potential with long-range hopping has re- 
ceived considerable attention recently. The Hamiltonian 
of such ID tight-binding model is expressed as 



H 



^e„ \n) (n\ + 

^ ^ J inn 

n n^m 



(9) 



where e n is energy level at nth site , uniformly distributed 



in interval [-W/2, W/2] and J mn = J/\m- 



(Jn 



0) is the long-range hopping amplitude. We will adopt J 
as energy units without loss of generality. The periodic 
boundary condition is applied. 

By using a supersymmctric method combined with a 
renormalization group analysis, Rodriguez et al. (2p| 
have shown the existence of extended states for energies 
within a range near the band top in one and two dimen- 
sional Anderson models. They found that MIT occurs 
only within the range of d < /i < 3d/2 in the thermo- 
dynamic limit, no matter how large the value of W is, 
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FIG. 3: Average concurrence iV (C 3 ) of the individual state in 
the long-range correlated disordered potential model as func- 
tions of energy for (a) a = 5.0, (b) a = 2.0 and (c) a = 1.0, 
respectively. The insets show concurrence for a typical dis- 
order realization at N=800. (d) Average concurrence N (C) 
over all states as functions of a. In all figures, from up to low, 
TV = 200, 400 and 800. 
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here d is the geometric dimensionality of the system. Us- 
ing finite size scaling analysis combined with the transfer 
matrix method, Xiong and Zhang [2l| found that there 
exists MIT at critical value \i c for some W. 

The concurrence N (C^^ as functions of energy are 
plotted in Fig. 4 (a), (b) and (c) for fi = 1.1, 1.5 and 
1.7 respectively. Here W — 5 is taken as example and 
the phenomena are similar for other W. The results 
are obtained for average of 200, 100, 50 samples [25[ for 
N = 200,400,800 respectively. When fi = 1.1 and 1.5, 
concurrence N (C 3 ) is large for states near the band top, 
which means there exist extended states. It is small for 
the states near the band bottom. This is quite different 
from that of above three models. For /i = 1.7, the con- 
currence is small for all states, which means all states are 
localized. 

The concurrence averaged over all states is plotted in 
Fig. 4(d). From the inset, we can get the inflexion for /i 
near 1.70, which is consistence with the upper limit for 
critical value /x c obtained in Ref.[2lJ. 

IV. CONCLUSION 

Using the measure of concurrence, mode entanglement 
sharing in one-particle states in four kinds of models 
is studied numerically. Concurrence N (C' 3 ) at a given 
state and N (C) averaged over all states are investigated. 
For (C' 5 ), the concurrence is large in extended states 
and small in localized states. There is a sharp transi- 
tion in the concurrence at mobility edge. N (C 3 ) gives 
the information about the localization behavior of the 
given eigenstate f3. From the curves of the N (C) vs 
the parameter of the models, which is A, V a — V&, a, 
and /i for slowly varying potential, random-dimer po- 
tential, long-range correlated disordered potential, and 
long-range hopping random potential, respectively, we 
can found clearly that there is an inflexion (or jump) at 
a critical parameter value, which is in accordance with 
that obtained by other methods. When parameter value 
is greater (smaller) than the critical parameter value, the 
system has only localized eigenstates, while when param- 
eter value is smaller (greater) than the critical parame- 
ter value, the system has both localized and delocalized 
states, which is different from that of one-dimensional 
Harper model. The inflexion or transition point in the 
curve of N (C) versus parameter of systems corresponds 
to the disappear of delocalized states. Therefore mode 
entanglement can be a new index to reflect MIT. 
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